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Abstract 

We propose a flexible framework for hedging a contingent claim by holding static positions 
in vanilla European calls, puts, bonds, and forwards. A model-free expression is derived for the 
optimal static hedging strategy that minimizes the expected squared hedging error subject to 
a cost constraint. The optimal hedge involves computing a number of expectations that reflect 
the dependence among the contingent claim and the hedging assets. We provide a general 
method for approximating these expectations analytically in a general Markov diffusion market. 
To illustrate the versatility of our approach, we present several numerical examples, including 
hedging path-dependent options and options written on a correlated asset. 
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1 Introduction 


Hedging derivatives using a static portfolio of standard financial instruments is a well-known al¬ 
ternative to dynamically hedging with the underlying asset. A static hedging portfolio is easy to 
construct and requires no continuous monitoring of the underlying or rebalancing over time. As 
such, static hedging strategies are more robust to significant underlying movements through mar¬ 
ket turbulence. Furthermore, static hedging portfolios are often useful fo r establishing no-arbitrage 
relati onships or bounds for exotic derivatives. This idea dates back to iBreeden and Litzenberger 
(119781) for standard options, and has been applied to exotic derivatives, such as basket options (see 


Hobson et al 


(12003)). 
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A fundamental result on static hedging due to lCarr and MadanI (|l998l l shows that any European- 
style claim written on a single underlying asset can be perfectly replicated by holding a fixed number 
of bonds and forwards, along with a basket of European calls and puts with the same underlying. 
The importance of this result is that it provides a model-free, perfect, static replicating strategy. As 
such, it also gives a no-arbitrage price relationship between the contingent claims and the hedging 
instruments. Nevertheless, there are also a number of limitations. In particular, the static hedging 
strategy requires an unbounded continuous strip of European calls and puts and must include the 
bond and forward in the portfolio. In reality, calls and puts are available only at discrete strikes in 
a finite interval. This leads to a practical question: how can one optimally construct a static hedge 
with only a finite number of calls and puts, with or without forwards on the same underlying? More 
generally, when there are simply not enough traded standard instruments t o achieve a perfec t static 
hedge, or when the hedger faces a binding cost constraint, the result of ICarr and MadanI (Il998l l 
provides no direction on how one might proceed. 

In this paper, we propose a flexible framework for hedging a contingent claim by choosing static 
positions in vanilla European calls, puts, bonds, and forwards. We are primarily interested in 
applications where the perfect static hedge is not available given a set of hedging instruments. To 
this end, we minimize the expected squared hedging error subject to a cost constraint. Our main 
result is a model-free expression for the optimal static hedging strategy, which involves computing 
a number of expectations that reflect the dependence among the contingent claim and the hedging 
assets. We provide a general method for approximating these expectations analytically in a general 
incomplete Markov diffusion setting that includes, but is not limited to, the well-known geometric 
Brownian motio n (GBM), Heston CEV a nd SABR models. 

Compared to lCarr and MadanI (119981 1. our framework includes a number of additional features. 
First, we allow for finite upper and lower bounds on the strikes of calls/puts used. Our static 
portfolio can involve any subset of the hedging assets among bonds, forwards, calls and puts, as 
opposed to include all of them. This gives the added flexibility to apply to underlying assets 
on which the forward contracts or some calls/puts are not written. Also, a cost constraint is 
incorporated into the hedging problem. When binding, this constraint may render a perfect static 
hedge impossible, and force the hedger to adjust the portfolio to minimize hedging error. While 
our methodology does not a priori assume the hedge is perfect, it can r ecover the perfect static 
hedge when it is available. This allows us to reconcile with the results in ICarr and MadanI (|l998l l 
as a special case of our fra mework. _ 

In the recent literature, I Carr and Wul (|2013l l work in a single-factor model and propose a finite 
approximation for the static hedging portfolio whose weights are computed based on the Gauss- 
Hermite quadrature rule. Also, there is a wealth of static hedging results specifically for barrier 
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options un d er on e- dimensional d i fFusion models; see 


Carr et al. 

(19981 

Carr and Lee 

(boog); 


Derman et al. 

(1995 


Carr anc 

Chou 

(1997) 


Carr and NadtochivI (j2nill ): 


Baxdosetal 


(j201(ll ). among 

others. In contrast to these works, our framework applies to other exotic derivatives and multi¬ 
dimensional diffusion models. We illustrate the static hedging strategies in three examples: Asian 
options, leveraged exchange-traded fund (LETF) options, and options with an illiquid underlying. 

The rest of this paper proceeds as follows: In Section [2l we formulate the optimal static hedging 
problem. In Section [3l we present our main results on hedging a contingent claim with a static 
portfolio of bonds, forwards and a strip of calls/puts. We also derive the optimal portfolio that 
consists of a finite set of assets. In both scenarios, we provide explicit, model-free optimal hedges. 
In Section 01 we discuss a practical method to numerically compute the hedging strategies in a 
general Markov diffusion setting. Lastly, in Section O we implement and illustrate our static 
hedging strategies in a number of applications. 


2 Problem formulation 

In the background, we fix a complete filtered probability space (12, T, (3“t)t>o, P), where P represents 
the physical probability measure and the filtration (3't)t>o represents the price history of the assets 
in the market. The market is assumed to be arbitrage-free but may be incomplete. We take as 
given an equivalent martingale (pricing) measure P ~ P, inferred from current market derivatives 
prices. For simplicity, we also assume a zero interest rate and no dividends. 

Our static hedging problem involves a group of hedging assets Z = {Z{x))x£i, with I being 
the index set. The number of hedging assets in I may be finite, countably infinite or uncountably 
infinite. The hedging assets could be, for example, bonds, stocks, calls, puts, forwards or other 
derivative securities. The price of each asset at any time t is denoted by Zt{x). 

We define a Static Portfolio as a signed measure 11 : 23(/) —^ P such that the static portfolio 
value at any time t is given by 

V^ = jn{dx)Zt{x). (2.1) 

In other words, n(dx) denotes the quantity of asset Z{x) of type x G dx held in the static portfolio. 
Observe that n(dx) may be negative, indicating a short position. Note that, while asset prices 
{Zt{x))x£i and the value of the static portfolio change with t, the number of units n(dx) 
remains constant for all t. 

Remark 2.1. We will consider two main examples in this manuscript: (i) hedging with calls/puts 
with strikes K in an interval K G [L, R), and (ii) hedging with a finite number of assets. In setting 
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(i), the signed measure II maps '?>{[L,R)) —?■ R. In this case, we will assume that 11 is absolutely 
continuous with respect to the Lebesgue measure and write n(ditr) = 'K[K)dK where vr is a function 
that maps {L,R) —)• R. In setting (ii), the signed measure II maps ®(Z+) —)• R. In this case, we 
will write n({i}) = Hi where vr. is a function that maps —)■ R. 

We now consider the contingent claim to be hedged at a future time T. Its market price at any 
time t is denoted by If the claim expires at time T, then H-p is the terminal payoff. We are 
primarily interested in situations where perfect static replication is impossible with a given set of 
hedging assets. Our goal is to minimize the expected squared hedging error of the static portfolio 
at time T subject to a possible cost constraint. We define the optimal static portfolio II* as the 
solution of the following optimization problem: 

n* :=argminE[(IljP-“r)% 5 := {R : Ijf < C}. (2.2) 

ne5 

That is, n* is the static portfolio that minimizes the expectation E[ {V^ — subject to the cost 
constraint < C. Note that the expectation in (|2.2p is evaluated under the physical probability 
measure P. Clearly, a perfect static hedge (V^* = Sy P-a.s.) is possible if and only if E[(I/^* — 
= 0. Note that the value of a portfolio at time t <T can be expressed as 

= jli{dx)nZT{x)\3^tl 

since all assets are martingales under the pricing measure P. Thus, the cost constraint Vq^ < C 
involves computation under the pricing measure P. 

Naturally, the optimal hedging performance and the corresponding static portfolio II* depend 
on the hedging assets available in the market, as well as the underlying price dynamics. Our main 
objective is twofold: (i) we provide a model-free expression for the optimal static hedging strategy 
when the hedging assets include bonds, forwards, vanilla European calls and puts on the same 
underlying; (ii) we discuss the implementation of the hedging strategies for a number of claims 
under Markovian diffusion dynamics. 

3 Methodology &: Main Results 

In this section, the set of hedging instruments contains a zero-coupon bond B, which pays one unit 
of currency at time T, a forward contract written on an underlying asset S with payoff {St — So), 
and T-maturity European puts and calls written on S. We assume there is a put at every strike 
K G [T,So) and call at every strike K G [So,R), with 0<L<So<R< oo. Let us denote by 
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g{K,ST) the payoff the call/put with strike K. That is 


g{K,ST) 


(K- 5 t)+ Ke[L,So) 
{St-K)+ KG[5o,i?) 


(3.1) 


el also 
(|l998ll 


While we observe L > 0 and i? < oo in practice and in our numerical e xamples, our inoc 
allows for L = 0 and i? = oo so that we can reconcile with the results in ICarr and Ma( 

(see Sect. 13.11 belowl. 

The terminal value of the static portfolio, composed of q bonds, p forwards and 'K{K)(iK units 
of European calls/puts with strikes in the interval dK, is given by 


cR 


Vf = q+p{ST-So)+ dK7riK)g{K,ST), 

where p, q and tt{K) may be either positive or negative (indicating a long or short position). 
The cost constraint is given by 


(3.2) 


H{TT,q):=q + J dK' k{K)z{K) < C, where z{K) :='E[g{K, St)]- (3.3) 

Note that, since the cost to enter a forward contract at inception is zero, the number of forward 
contracts p in the static portfolio plays no role in the cost constraint. 

With given by (|3.2I) . and cost constraint by (13.3p . algebraic calculations show that the static 
hedging problem (12.2p is equivalent to solving for 


where 


{7T*,q*,p*) := argmin J(7r,g,p), 
(w,q,p)eS 


S := {{Tr,q,p) : H{7r,q) < C}. (3.4) 


J{7T,q,p) :=E[{Vf-E t)^] 

pR pR 

= q‘^+p‘^j:+ / dKdK'7r{K)'ip{K,K')Tr{K') 

J L J L 

rR 

+ 2qpl3 + 2q dK 7r(K)z(K) — 2q^ 
pR pR 

+ 2p dK T{K)y{K) - 2pe - 2 I dKT{K)-i{K), 

J L J L 

and we have defined the expectations: 

/3 = E[5r-5o], e = E[iST-So)ET], E = E[iST-So)% 

e = E[HT], ^{K) = E[ETg{K,ST)], ij{K,K') = E[g{K,ST)g{K',S t)], > 

z{K) = E [g{K, ,St)], y{K) = E [{St - So)giK, 5t)]. 
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(3.5) 


In order to state and prove the optimal hedging strategy in this setting, we need the following 
Lemma. As preparation, it is convenient to introduce the probability density functions of St under 
the physical (i.e., historical) and risk-neutral probability measures 

rs^{K)dK = P{St edK), fs^{K)dK = P{STedK). (3.6) 


Lemma 3.1. Assume the random variable St has a strictly positive density Tsj, & (7^(11+). Recall 
the function if as defined in (133]), and let f : R+ —>• R 6e C'^(R+). Then the solution vr of the 
integral equation 


is given by 


f{K) = j\K'7T{K')ifiK,Kfi, 


<K) := dl 


\TsAK))' 


(3.7) 


(3.8) 


Proof. In what follows, let 11 be the anti-derivative of vr and 11 be the anti-derivative of 11 so that 
n' = TT and = vr. Observe from ()3.1|1 that 


0= lim g(K'K)= lim dK'g{K'-,K)= lim g(K'K)= lim dK'g{K',K)- 

K'^L ^ ' K'^L ^ ' K'^R ^ ' K'^R ^ ' 

Let us further observe that d\g{K, s) = 5{K — s). Then, equation (|3.7p implies 


dlf{K) = dl p dK' T{K')if{K,K') 

rR 


= ^ dK' 7r(iL')5|:E [g{K, ST)g{K', 5^)] 

pR roo 

= / dK'Ti{K')dl / dsg{K,s)g{K',s)Ts.r{s) 
Jl Jo 

rR roo 

= / dK'7r{K') / dsdlg{K,s)giK',s)TsM 

Jl Jo 

rR roo 

= / dA:'7r(A:') / ds<5(A:-s)c/(A:',s)rs^(s) 

Jl Jo 

rR 

= rsAK) dK'T{K')g{K',K) 

= Ts^ {K) (g{K', K)U{K') - dK'g{K', A)n( A') 


rR 


+ 


dK'li{K')dl,g{K',K) 


(by (I33|)) 


(3.9) 


(integrate by parts) 
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= rsAK)niK) 


(by (1331)). 


(3.10) 


To obtain ()3.8p . simply divide (13.101) by Tsj^iK), differentiate both sides twice, and use 11 = vr. □ 

Using Lemma EH we can now state and prove the optimal hedging strategy. To this end, we 
define the function 


7r{K, A) := 


dl^jK) X TsAK) 
TsAK) 2TsAK) 


(3.11) 


where K G [L,R], A G M, the densities Tsj, and Tg^, are defined in ()3.6I) and the function 'y{K) is 
given in ()3.5p . 

Theorem 3.2. Assume the random variable St has a strictly positive density Tsj, G C'^(R+) under 
P and a density TSj, G C'^(R+) under P. Assume further that 7 G C'^(R+). Finally, assume the 
matrix inverses defined in ()3.12p and (|3.13p are well defined. Then the optimal strategy {TT*,q*,p*) 
that solves the optimal static hedging problem dH is given by 


A*,q*,P*)={ 


(7r(.,A^),g^,p^) if q^+l dKn{K, XAz{K) < C 


r-R 




{Tri-,X^),q^,p^) else, 


where 


X^ := 0 , 




and 






V 


\P 


1 /3 


U 


'1 /3y7«-U'lA'^(A')8|-(g^) 

7 U-/«dA'!,(A)aS, 


(3,12) 


2 (Tst{K)\\ 


-1 


1 0 


1 fR 

2 


^SriK) 
'K\ TFAK) 


1 rR 


J^dKAK)dl 


(^-lL<^KziK)d],(p0^)\ 


0 - dK yiK)dl{p0^) 

[c-lLdKAK)dl{Pl^)) 


(3.13) 


Proof. First, we define the Lagrangian associated with (|3.4p : 


LA,Q,pA) ■= JA,q,p) - A {H{TT,q) - C) 

= q^ + — 2p6 — 2qf^ + 2qpfi — Xq + XC 

pR pR 

+ / / dKdK’Tr{K)fi{K,K')Tr{K') 

Jl J l 
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+ f dK7r{K)(2qz{K)+2py{K)-2-f{K)-Xz 


where L{-,q,p,\) acts on functions in C{[L,R]). The Karush-Kuhn-Tucker (KKT) conditions, 
necessary for optimality, are (below, rj is an arbitrary C{[L,R]) function satisfying \\r]\\^ < oo) 


stationarity: 


d 


0 = —L{7T + er],q,X)\^^^ 


= 2 


€=0 


rR 


dKp{K)(^l^qz{K)+py{K)--/{K)-^^K)^ + j dK'7r{K')'iP{K, K' 


rR 


0 = qz{K) + py{K) - 7 (K) - ^z{K) + J dK' tt{K')xIj{K, K'), 
stationarity: 0 = dqL{TT, q,p, X) 

rR 

= 2q — 2^ + 2pf3 — X + 2 J dK 7r{x)z{K), 
stationarity: 0 = dpL{TT, q,p, A) 

rR 

= 2pX:-29 + 2qp + 2 J dK7r{x)y{K), 
comp, slackness: 0 = A • q) — C) 

= X-(q+[ dK^T{K)^K)-C 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


Note that (j3.14D is of the form (j3.7l) with f{K) = ^{K) + ^z{K) — qz{K) —py{K). Thus, using 
Lemma l3.ll we obtain 

7 r(iL) = d'i 

Next, noticing that 


,2 , dl^{K) + ^dlz{K) - qdlz{K) - pdly{K) 


Tst{K) 


(3.18) 


^ 1*00 ^ 
dlz{K) = alE [g{K, 5r)] ds g{K, s)Tsr (s) 

Jo 

/*oo ^ /•oo ^ ^ 

= / dsdlgiK,s)TsM= ds5{s-K)TsAs) = rsAK), 

Jo Jo 

poo 

dlz{K) = [g{K, 5r)] =8], ds g{K, s)Ts^ (s) 

Jo 

poo poo 

= dsd],giK,s)Ts^{s)= ds5{s-K)TsAs) = rsAK), 

Jo Jo 

poo 

dly{K) = 5iE [{St - So)g{K, 5t)] = / ds (s - So)g{K, s)rs^(s) 

poo poo 

= ds{s - So)dKg{K,s)TsT{s) = ds {s - So)S{s - K)rsT{s) 
Jo Jo 





{K-So)rsAK) 


and substituting these expressions into (13.181) . we see that Tr{K) in ()3.18p coincided with the ex¬ 
pression given in (j3.1ip . Next, inserting expression (j3.1ip into the KKT conditions, (j3.15p . ()3.16p . 
and (I3.17p . gives the following system of three equations 


0 = 2q 


fR 

2C + 2p/3-X + 2 J dKz{K)d‘j^ 


[ ^sAK) ) ’ 


cR 

0 = 2pS - 20 + 2g/3 + 2 J dK y{K)d‘i 


f dlAK) + k^sAK) 
\ ^sAK) 


0 = A • 



r-K 


dKAK)dl 


( d],AK) + ^rsAK) 

V ^sAK) 



(3.19) 

(3.20) 


The above system has two possible solutions corresponding respectively to the cases A = 0 and 
A / 0. For the case A = 0, the triplet (j),q,X) is given by (j3.12p . On the other hand, when 
A / 0, the triplet {p,q,X) is given by (13.131) . Finally, the KKT conditions are necessary conditions. 
Since both the objective function and constraint are convex, and the primal pr oblem is feasible 
(Slater’s condition), the KKT conditions are also sufficient for optimality (see, e.g. IZalinescul (|2002l . 
Theorem 2.9.3)). □ 


In Theorem 13.21 the solution A*,q*,p*) = X^),q^,p^) corresponds to the unconstrained 

optimization problem. But if the associated cost is less than C, that is, 

pR 

+ j dK7r(K, X^)z{K) < C, 

then (7r(, •, A^), q^ ,p^) must coincide with the solution X'^),q‘^,p^) of the constrained opti¬ 
mization problem with initial cost (upper bound) C. On the other hand, if the unconstrained 
optimization problem admits a cost greater than C, then the corresponding constrained optimiza¬ 
tion problem has the solution (7r(-, A*^), g^,p^) where, by construction the constraint is binding: 
q^ + dKAK, X^)z{K) = C. 

In fact, from (13.lip , the first term in parenthesis in the optimal strategy vr* can be interpreted 
as a conditional expectation. Heuristically, we have 

10A " fe'(Sr)2Tl = E[=t|St = K], (3.21) 

In other words, the number of units of call/put held at strike K involves computing the expected 
terminal claim Et conditioned on the terminal price of the underlying taking value K. 
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Remark 3.3. Although we considered hedging with European call/puts on a single asset S, the 
results of this section can be extended to the case where one hedges with European calls/puts on 
n assets , S'^,... S'^. The only difficulty that may arise is in solving the equations that result by 
imposing the KKT conditions. 


3.1 Connection to 


Carr and MadanI ( 19981 ) 


Let us recall the main result in 


Carr and Madan (1998): if / : 1R+ —>• R satisfies / G C^(R+), then 


p(X> 

fiSr) = f{So) + f'iSo){ST -So)+ dKr{K)g{K,ST). (3.22) 

Jo 

As such, a contingent claim with payoff /{St) can be perfectly hedged by holding /{Sq) bonds, 
/'(So) forward contracts and a basket of puts and calls, where the weight of the put/call with 
strike K is f"{K). The following corollary proves that equation (|3.22l) is indeed a special case of 
our Theorem 13.21 

Corollary 3.4. Consider a European-style contingent claim with payoff Et = /{St) as in (I3.22p . 
Assume f G C^(R+). Let L = 0 and R = oo so that a call/put on S is available at every strike 
K G [0, oo). Then, under the assumptions of Theorem \3.^ the optimal static hedging portfolio with 
no cost constraint is given by ('/r*,g*,p*) = {Tr{-, X^),q^,p^) where 

7r{K, A^) = f'{K), = 0, = /(Sq), p^ = /'(Sq). (3.23) 


Proof. We must show that {Tr{-, X^),q^,p^), given by (|3.23|) . satisfies (|3.19l) and (l3.2Up . Eirst, we 
observe that 

/•OO 

dl-f{K) = dlE[g{K,ST)ET] = d\ / ds g{K, s) f {s)VsM 

/•oo /•oo 

= / dsdlg{K,.s)f{.s)TsAs)= ds5{s-K)f{.s)TsAs) = f{K)rsAK). 

Jo Jo 

Thus, using A^ = 0, we see from (|3.1ip that 


T{K,\^) = d'i, 


l{K) \ 
^sAK)) 


rP 


- rP 


^5,(iL)/(iL)^ _ 




= r{K). 


Next, dividing equation (I3.19p by two and rearranging terms we find 


E/(Sr) = g + pE(5T-So)+ / dKE[g{K, ST)]f{K). 

Jo 

This equation will clearly be satisfied if q = /(Sq) and p = /'(Sq). To see this, simply take the 
expectation of (I3.22p . Next, dividing equation (|3.2Up by two, and rearranging terms we have 


E [{St - So)/{S t)] = qE {St - So) + pE [(5t - 5o)"] + / diLE[(5T - So)g{K, St)]/'/^). 
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This equation will also be satisfied if g = /{Sq) and p = f'{So). To see this, simply multiply (j3.22l) 
by {St — Sq) and take an expectation. □ 


3.2 Connection to static hedging with finite assets 

Our framework can be related to static hedging with a finite number of assets. In this case, the set 
I has a finite number N units of hedging assets, and the static portfolio value can be expressed 
as a finite sum: 

N 

= Y,^iZT{i). (3.24) 

This is indeed the discrete version of the static portfolio in (12.ip . 

Assumption 3.5. We assume that the random variables {ZT{i))i£i are elements of L^(P) and are 
linearly independent. Stated in financial terms, this assump tion simply requires that none of the 


hedging instruments is redundant, as defined in (jPuffiel . l200ll . Chaper 2). 


With given by (I3.24p . a direct computation shows that the static hedging problem amounts 
to determining the optimal strategy 


TT* := argmin J(7r), 
ttGS 


S :={tt: H{tt) < C], 


(3.25) 


where J{'k) and H{'k) are given by 

•^( tt ) = ^ ^ TTiA,jT^j - 2 ^ 


H{7r) = 




and 


= IE [ZT{i)ZT{j)], 


7i = E [ZT{i) Ht], 


Zi = E [ZT{i)]. 


(3.26) 


Compared to the “continuous” case in (I3.2p - (l3.5p . the objective function again involves the 
expectations of products of payoffs, namely, E [ZT{i) and E [Zxii) ^t]- Note that in this discrete 
case we can consider general claims, not limited to forwards, puts, and calls, and continue to derive 
explicitly the optimal static portfolio. 

Proposition 3.6. Let V’ be the square matrix whose {i,j)-th component is Then, under 

Assumvtion \3.5l the optimal static portfolio vr*, defined in (|3.25l) . is given by 


* 

7T = < 


vr^ = 


if z'^TT^ < C, 
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else. 


(3.27) 













where ^ is the inverse of if:, and z, vr and 7 are column vectors whose i-th components are zi, TTj 
and 77 respectively. 

We provide a proof in Appendix |Al The vector corresponds to the optimal strategy without 
a cost constraint. If the unconstrained optimization problem has a cost < C, then is the 

optimal strategy for the constrained optimization problem. On the other hand, if the unconstrained 
optimization problem has a cost > C, then the solution of the constrained optimization 

problem is given by which, by construction, has a cost equal to C, that is, = C. Lastly, 

we emphasize that the optimal static hedging strategy with discrete strikes can be quite different 
than the optimal strategy when continuous strikes available but implemented at discretized strikes. 
We will visualize the difference in Section O 


Remark 3 . 7 (R elation to Markowitz mean-variance portfolio optimization). In his seminal work, 


Markowita (jl952l l solves the problem of minimizing portfolio variance for a given level of expected 


return. Mathematically, the minimization problem is given by 


min w^TiW, W := {w : pL^w > m and ||t(;|| = 1}, (3.28) 

where w are the portfolio weights to be found, S and [i are, respectively, the covariance matrix and 
expected returns of a group of assets, and m is the minimum level of expected return. Interestingly, 
the portfolio optimization problem (|3.28p has the same structure as the static hedging problem 
()3.25p . which, in matrix notation, is given by 

min , 5 := {tt : T'^tt < C}. (3.29) 

7rG<S 

Though, clearly, the economic interpretations of ()3.28p and ()3.29l) are distinct. 


4 Implementation Under a Markov Diffusion Framework 

Thus far, we have made no assumption about the dynamics of the underlying S. To illustrate 
the performance of our static hedging strategies, we now present the calculations and numerical 
implementation under a general incomplete Markov diffusion market. The analytic approximations 
we present below are useful when the claim to be hedged is European-style. Specifically, the 
payoff r.T may be some function h of the final value of a d-dimensional Markov diffusion X. Note, 
by allowing components of X to be the quadratic variation or running average of other components, 
our definition of European-style claims allows for path dependence and includes both Asian options 
and options on variance/volatility (e.g., a variance swap). Extending the approximations to cases 
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where H-r is a barrier-style claim or look-back option is not trivial and is well beyond the scope of 
this paper. 

Let X = (X^,..., G IR'^ be a Markov diffusion satisfying the following stochastic 
differential equations (SDEs) under P and P, respectively 


dXt = n{t,Xt)dt + a{t,Xt)dWt, (under P) (4.1) 

dXt = Jl{t,Xt)dt + a{t,Xt)dWt. (under P) (4.2) 

Here, W (resp. W) is an m-dimensional Brownian motion under P (resp. P), and the functions 
J1 and a map 

/i : P+ X J1:R+ a : E+ x 


Let us suppose that the terminal values of the hedging assets {Z{i))i^i from Section 
S from Section [3] and the claim H to be hedged are given by 


the stock 


gi{X)r) = g{Ki,e^T)^ 


St = \ogX]p, 


It = KXt), 


where the function h maps E'^ i—>• E and for each i the function gi maps E E+. Since S = logX^ 
is traded, in order to preclude arbitrage, we must have 


= 


m 


In order to implement the optimal hedging strategies (Theorems 13.21 and 13.6p . we must com¬ 
pute the expectations defined in equations (13.51) and (I3.26p . For general dynamics of the form 
(jlll-dlSI), closed-form expressions for these expectations are not available. Moreover, comput¬ 
ing these expectations via Monte Carlo simulation is not practical, since, in the case of Theorem 
Ea the expectations appear in the integrands of various integrals. As such, we provide here a 
method for obtaining analytic approximations the expectations in (13.51) a nd (13.261). The methods 
that w e describe below were developed first formally in a scalar setting in 


Pagliarani and 


(I 2 OI 2 II and later extende d to multiple dimensions with rigorous error bounds in lLorig et al 


and 


ascucci 


(I2ni5bl l 

Lorig et al.l (j2015al L Here, we give a concise review of these methods and also provide some 


extensions, which are needed to implement Theorems 13.21 and 13.61 

We fix a time T > 0 and consider an expectation of the general form 


u{t,x) = E[(^(Xr)|Xt = x]. 


t < T. 


(4.3) 


Under mild conditions on the drift /r, diffusion coefficient a and terminal data ip the function u 
satisfies the Kolmogorov backward equation. Omitting x-dependence below to ease notation, we 
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have 


{dt +A{t))u{t) = 0, 


u{T) = if, 


(4.4) 


where A{t) is the generator of X under probability measure P. Explicitly, the operator A{t) is 
given by 

d d d 

yi(t) = =: ^ a„(t,x)5“, (4.5) 


2=1 2=1 j = l 

where we have introduced standard multi-index notation 


l<|o|<2 


a — («!, • • • , Ud) G INq, 


a = 




d^ = tld^. 


2=1 


2=1 


Remark 4.1. To compute u{t,x) := Wi[ip(XT)\Xt = x], one would simply replace A{t) in (14.41) 
with 

d d d 

A{t) = + =■ a„(t,x)5“. 

i=l i=l j=l l<|a|<2 

Our goal is to hnd an approximate solution to PDE (14.4j) . thereby obtaining an approximation for 
the expectation ()4.3I) . To this end, we expand each coefficient as a Taylor series about a fixed 
point X G R'^: 


OO 

G'a — ^ ^ 0^CK,n; 
n=0 


, S dxaa{t,x) _ g 

y,n{t,x) := 2^ - p - (x-x)^ 


x^ = xf*. 
2=1 


(4.6) 


Here, we assume implicitly that the coefficients Oq, are analytic. However, we will see in Definition 
I4.2l that the A^th-order approximation of u requires only that the coefficients be C'^(R‘^). Combining 
()4.5p with (14.6p we see that the operator A{t) can be written as 


A{t) = Ao{t)+ ‘Bi{t), = y^^Anjt), 


22=1 


Inserting (14.7p into PDE (14.411 we have 

{dt +Ao{t))u{t) = -‘Bi{t)u{t), 

and, hence, by Duhamel’s principle 


‘^n(t) = Y «a,n(i,a;)9“. (4.7) 

l<|o|<2 


u{T) = ip, 


u 


(t) =‘Po{t,T)ip + / dti‘P{t,ti)'Bi{ti)u{ti), 


(4.8) 
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where J’o(i, T) is the semigroup generated by -do(t). Explicitly, we have 


•Po{t,T)tp{x)= dy To{t,x-,T,y)ip{y), 


( 4 . 9 ) 


where To(t, x',T, ■) is a Gaussian kernel whose mean vector m(t, T) and covariance matrix C{t,T) 


are 


in{t,T):—x + J^ ds ^ 0 ( 1 ^ 0 ,••• ,o),o('S) ,o),o('S) 

^2a(2,o,--- ,o),o('S) “(1,1,••• ,o),o('S) 

“(i,i,--- ,o),o('S) 2 a(o, 2 ,--- ,o),o('5) 


C(t,r) :=^ 


ds 


\ “(i,o,--- ,i),o('S) “( 0 , 1 ,■■■ ,i),o('S) 


“(o,o,-,i),o(s)) > 

“( 0 , 0 ,••• ,l),o('S) \ 

“( 0 , 1 ,,i),o('S) 


2 a(o,o,---,2),o('5)/ 


Observing that u appears on both the left and right-hand side of (14.81) . we iterate this expression 
to obtain 


u{t) = 7o{t,T)^p+ / dti / dt 2 ' 


dtk 


k=i • 


'tk-l 


'?o{to,ti)'Bi{ti)‘J’o(ti,t2)‘Bi{t2) ■ ■ ■ '?o{tk-i,tk)'Bi{tk)7o{tk,T)ip 


= J’o(i 


oo n pT 

EE/ dti / dt2 ■ 

n=l k=l •^*1 


dtk 


'tk-l 


'?o{t,ti)Aij^{ti)‘?o{ti,t2)Ai2{t2) • • • ?o{tk-i,tk)Ai^{tk)'?o{tk,T)^, (4.10) 

In,k = {* = (ii, *2, • • • , 4) G : n -b ^2 H-b 4 = n}. (4.11) 

where, in the second equality, we have used the fact that (Bi(t) is an infinite sum and we have 
partitioned on the sum of the subscripts of the Ai{t) operators. Expression (|4.10l) motivates the 
following definition. 

Definition 4.2. Let u be the unique classical solution of (14.41) . Assume the coefficients aa{t, •) are 
C'^(1R'^) for all t G [0,T]. Then, the Nth order approximation of u, denoted by un, is defined as 


N 

Un •— ^ ^ llnj 
n=0 


Uoit) := yo{t,T)(p, 


where '?o{t,T) is the semigroup generated by Ao(t) and 

n ,.T .T i-T 

Un{t):='^ dti dt2--- dtk 
Jt J t\ J 
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‘yo{t,ti)Ai^{ti)7Q{ti,t2)Ai^{t2) ■ ■ ■ CPo(4-i,4)-Aifc(4)J’o(4)(4-12) 

'^^^71,k 

The Nth order approximation of the transition density TM{t,x-,T,y) = '^n=o^'ri{t,x;T,y) is ob¬ 
tained by setting the terminal data equal to a d-dimensional Dirac mass (p = 5y. 


Recall from (I4.9p that the semigroup operators ‘?o{ti,tj) are integral operators. Thus, as written, 
expression ()4.12l) is difficult to evaluate. The following Theorem shows the nth order term Un can 
be easily computed as a differential operator Ln acting on uq. 


Theorem 4.3. iLoria et al\.\2015d. Thorem 2.6) Let {un)n>o be as given in Definition (I4.2p . Then 


Unit) = Lnit,T)uoit), 


n T pT 

Lnit,T) = '^ / dfi / dt2--- 

k=l 



where In,k is defined in (14.1111 and 


9iit,tj):= 

l<|a|<2 


X(t, tj) := X -|- m(t, tj) + C(t, tj)Vx. 


Theorem 14.31 provides provides a method of approximating the expectations in (13.51) and (I3.26P 
analytically. The expectations in (13.261) are sufficient to implement Theorem l3.6l which can be used 
compute optimal static hedges with a discrete number of hedging assets. 

In order to compute optimal hedges in the case with a strip of calls/puts through Theorem 13.21 
in addition to computing expectations in (|3.5p . we need numerically evaluate 


^sAK) 

^sAK) 


and 


dhiK) 
^sAK) ’ 


where Tsj, and Tsj, are the densities of St under P and P, respectively, and 7 is defined in (j3.5p . 
To this end, we compute 

fsAK) f^i(log/^) 

TsAk) K-^r^,^iiogK) r;,i(iogi^)’ 


where T^^^ and 

dlAK) 

^sAK) 


1 are the densities of X)p = logS^ under P and P, respectively. We also have 

,^-^-E[dlg{K,e^T)h{XT)] 


^sAK) 


^[5K{e^T)h{XT)] 


(logR) 


dy 5KA^'')hiy)TxAv) 
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= [ dy' h{{log K,y')) ^ , y'= {y 2 ,y 3 ,... ,yd), (4.13) 

where Fxy is the density of Xt under P. Note that (j4.13D is in fact Wi[h{X'j’)\X^ = logiF] (see 

(I32I])). 


4.1 Accuracy results 

In order to establish the accuracy of approximation given in Definition 14.21 we must make a few 
assumptions about the coefficients (uq,). 

Assumption 4.4. There exists a positive constant M such that the following holds: 

1. Uniform ellipticity: < Yl,a= 2 ^a{i,x)^°^ < M|^p, Vt G P+, G P'^, 

2. Regularity and boundedness: the coefficients aa{t,-) G C'^"'“^(P‘^) and partial derivatives 
d^Oait, •) with |/3| < N are bounded by M for all t G P+. 

We begin by establishing the accuracy of the approximation uw 

Theorem 4.5. Let Assumption \4-4\ hold and fix x = x. Suppose the terminal datum ip is at most 
exponentially growing and p G C'^“^(P'^) for some 0 < k < 2, where C“^(P'^) denotes the space of 
functions that are not necessarily continuous. Then we have 

AT+l + fc 

\u{t, x) — U]\[{t,x)\ < C{T — t) 2 j 


where C is a positive constant the depends on M, N, the terminal datum p and x. 


Proof. See (|Lorig et all . l2015al . Theorem 3.10 and Remark 3.11). 


□ 


Next, we will derive an approximation for df;^Tsj.{K), which is needed to compute (13.111) . and 
establish the accuracy of this approximation. To this end, let us define 


cit,x;T,K) = E[{ST-K)+\Xt = x], 
Since X^ = log S we clearly have 


u{t, x; T, k) := E [(e^x — e*^)’''| = x] 


c{t, x; T, K) = u{t, x; T, k{K)), k{K) := log A. 

Note that the density of St can be obtained by differentiating c(t, x; T, K) twice with respect to 
strike 


F{St G dK\Xt = x) = Tst {t, x; T, K)dK = d]^c{t, x; T, K)dK 
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= d^Ku{t,x;T,k{K))dK. 
Thus, we define Tsj,,n, the Nth order approximation ofTsj., as 


(4.14) 


rsj,,Ar(t,x;r,iir) ■.= d\uN{t,x\T,k{K)). (4.15) 

The following theorem gives the accnracy of d’^TS t AK) 

Theorem 4.6. Let Assumption \4-4\ hold and fix x = x. Then we have 

\d^rsAt,x-,T,K)-d^rs^,N{Lx-,T,K)\<C{T-tA-^y\ (4.16) 


where C is a positive constant the depends on M, N, K and x. 
Proof. Prom (j4.14p and (|4.15|) . we have 


df, {t, x; T, K) - Ts^,N{t, x; T, iP)) 

.. n+2 

i=i 


where the coefficients (aA'^) a r e inte gers whose precise value is not important. Next, by Theorem 


4.4 of 


Pagliarani and Pascuccil (|2ni5l ). we have 


\d]iu{t,x-,T,k) - d^UN{t,x-,T,k)\ < C(r-t)(^-”+2)/2, 


(4.18) 


where the constant C depends on M, N, k and x. The accuracy result (j4.16p follows directly from 
(HTTP and (liTHI) . □ 

A similar analysis yields the same order of accuracy for the approximations of dffiTS tA) 
df^'y{K). For brevity, we do not repeat the computations here. 


5 Examples 

In this section, we compute the optimal static hedges, along with sensitivity analysis, in a variety of 
practical applications. The resulting static portfolio payoff profiles are also shown to demonstrate 
the effectiveness of our methodology. 
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5.1 Hedging an option with options on a correlated underlying 

In our first example, we consider two correlated underlyings S = and V = e^^, where the 
log-price pair X = satisfies 


dXl = ^il{t,Xt)dt + ai{t,Xt)dW^\ 
dX^ = Mt,Xt)dt + a2{t,Xt){pdWl + 

/ 

dX^ = Mt,Xt)dt + ai{t,Xt)dWl, 

dX? = Mt, Xt)dt + a2{t, Xt){pdWl + ^l-p^Wt^).^ 


(under P) (5-1) 


(under P) 


We assume that S is traded and is therefore a martingale under P. As such, we must have 
Jli = —^crf. Although we do not assume it, if V is traded, it too must be a martingale under 
P and in this case the drift ^2 must satisfy JI 2 = — 

Suppose now we have sold a European call written on V with maturity date T and strike price 
K'. This is the claim to hedge, so we denote 


Et = {Vt - K')+ = (e^T - K')+ =: h(K',X^). 


(5.2) 


We wish to statically hedge this option with bonds, forward contracts and European calls/puts 
written not on V but on S. This situation may arise for a variety of reasons. First, it could be that 
neither V nor options on V are liquidly traded, or the options are traded at very few strikes, as 
is common in the commodity markets. Second, even if V is traded, an investor may be prohibited 
from trading V and options written on V for legal reasons. Third, the market for options on S may 
simply have superior liquidity as compared to V. For instance, one may resort to forwards and/or 
options on the S&PSOO index in order to hedge an option written on a component or non-component 
stock. 


Continuum of Strikes K G [L, R) 

Let us recall the setting of Section[21 in which the bond, forward on S, and calls/puts are available 
at every strike K G [L,R) (see ()3.ip L To compute the optimal hedging strategy (vr*,g*,p*) using 
Theorem 13.21 we require expressions for the expectations in (13.51) . For the current application, they 
are given by 

S = E[(e^T_e")2], i;{K,K') = E[g{K,e’^T)9{K',e’^T)]^ /? = E [(e^r _ e")], 

z{K)=M[giK,e^h], yiK) = E[ie^T - e-)g{K,e^h], ^ = E[h{K',X^)], 

e = E[{^^T-e-)h{K\Xl)l j{K) = E[g{K,e^hh{K',X^)], z{K) = E[g{K,e^hl 
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where g is defined in ()3.1I) and h is defined in (I5.2p . For certain model dynamics, these expectations 
can be computed explicitly. In cases in which the expectations cannot be explicitly computed, 
analytic approximations can be obtained using Theorem 14.31 

In Figured] we consider the case where and are correlated arithmetic Brownian motions. 
We assume that both S and V are traded assets, which requires that Jli = and ^2 = —^(^ 2 - 

We illustrate the optimal static portfolio {q*,p* and examine the effects of the correlation p 
and the cost constraint C. In addition to plotting the units of calls/puts held, that is, 'k*{K) as a 
function of strike K, we also plot the static portfolio’s terminal value 

[R 

d>(5r) := q*+P*{ST - So) + dKTT*{K)g{K,ST) (5.3) 

as a function St, and compare this portfolio profile against the payoff {Vt — K')'^ . 

From Figured) we see two clear effects as the correlation increases to I. First, the density of 
calls/puts in the optimal hedging portfolio becomes more concentrated near K = K' . Second, the 
payoff function <I> more closely matches the payoff function h{K',log{-)) = (• — K')^ of the option 
to be hedged. This is intuitively what one would expect since, if S and V are perfectly correlated 
(i.e., p = 1 ), then the two assets are identical, i.e. St = Vt, so holding the call on St with strike K' 
will be the optimal and perfect hedge. Somewhat less intuitive is the effect of the cost constraint 
on the optimal hedging portfolio vr*. In Figured) we see that, as the cost constraint becomes more 
severe (i.e., as C decreases), the optimal static hedging strategy is to sell more puts and calls with 
strikes K far away from K' . Intuitively, selling those options helps reduce the hedging cost in order 
to satisfy the cost constraint while not increasing the hedging errors around K'. 

Hedging with discrete strikes (Ki) 

As in Section [321 we now assume that European-style calls/puts are available only at a discrete 
strikes {Ki)i^T We set Zt{0) = Bt = 1 and Zt{i) = W,[g{Ki, ST)\3't] for i ^ 0. In order to compute 
the optimal static portfolio vr* using Theorem 13.61 we compute the expectations in (|3.26l) . which 
are given by 

" 00,0 = 1 , i>o,i = i>i,o = '^[g{Ki,e^T)] = E[g{Ki,e^T)g(^Kj,e^T)]^ 

70 = E[h{K',X^)], = E[g{K„e^T)h{K',X^)], zo = li = 1 , 

Zi = Elg(Ki,e^T)j^ = E lg(Ki, iJ^O. 

For certain model dynamics, the above expectations can be computed explicitly. In cases in which 
the above expectations cannot be computed explicitly, they can be approximated analytically in a 
Markovian framework using Theorem 14.31 
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As in the previous case with continuum of strikes, we assume that the log prices and X'^ 
are correlated arithmetic Brownian motions. In Figure [2] (bottom), we show the optimal static 
portfolio weights and illustrate the effect of the correlation parameter p and cost constraint C. In 
addition, we plot the static portfolio profile (terminal value as a function of St), denoted by 

^St) := TTo* + Y,<9{K^,ST). (5.4) 

i 

As the correlation increases from 0.5 to 0.9, we see that the portfolio profile matches more closely 
the claim payoff h{K' ,\og{-)) = (• — K')'^, especially for large and small values of St- However, in 
contrast to the case with continuous strikes, a perfect match is not possible due to the availability 
of only a limited number of strikes. On the other hand, the effect of the cost constraint on vr is 
less straightforward. Specifically, with a more stringent cost constraint (i.e., C —>■ 0), the optimal 
static portfolio tends to have a negative payoff for low strikes and a positive and increasing payoff 
for high strikes, resulting in a higher quadratic hedging error. Furthermore, the resulting portfolio 
profile is lower for all realization of St when the allowed portfolio cost C is reduced. 

Moreover, under all conditions we observe that the magnitude of tt is greatest for strikes Ki 
near K'. However, due to the fact that we have a discrete number of hedging assets, the sign of 
vr oscillates as a function of strike. In contrast, when using the optimal static strategy derived 
from the model with a continuous strip of calls/puts, the sign of the optimal density vr does not 
oscillate. This has important practical consequences. From the hedger’s point of view, it is far 
easier to take only or mostly long positions, rather than taking alternating long and short positions 
in options. Therefore, even though there is a finite number of strikes in practice, one can also adopt 
a discretized version of the optimal continuous density in Figure [T] (top), as an alternative to the 
optimal discrete strategy vr(iFj) in Figure [2] (top). Hence, by discretizing the optimal static hedging 
strategy in Theorem 13.21 (which assumes options trade with strikes in a continuum) one obtains a 
useful alternative to the strategy derived from assuming discrete strikes from the onset. 

5.2 Hedging an LETF option with options on the reference 

An exchange traded fund (ETF) is typically designed to track a reference index, denoted by S' = e^^. 
The reference index dynamics under under the physical measure P and pricing measure P are given 
by two-dimensional SDE (|5.1I) . where the second component V = log now represents the driver 
of volatility rather than a traded asset. Now, we introduce L, a leveraged exchange traded fund 
(LETE). An LETE is a managed portfolio that returns a pre-specified multiple i of the daily return 
of the reference index S. To illustrate this in the simplest form, with zero management fee, interest 
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and dividend rates, the dynamics of an LETF with leverage ratio I are related S as follows: 


Lt St ■ 


The m ost common values for the leverage ratio i are {—3, —2, —1, 2,3}. As lAvellaneda and Zhang 
(120101 ) show, the value of the LETF L at any time T > 0 is given by 


Lt fSTV 


^ crlit,Xt)dt^ 
= exp - Xl) + . 


(5.5) 


where {X^)t is the quadratic variation of X^ over the interval [0,T]. Observe that the value of Lt 
depends not only on the terminal value St = e^r^ but also on the integrated variance (quadratic 
variation) of log S = X^. Thus, the value of L depends on the entire path of S. 

Options on LETFs of different leverage ratios are also widely traded on the Chicago Board 


Options Exchange (CBO E). The pric i ng of t 
been recently studied in iLeung et ahl (j2014l i: 


lese options unde r stoch astic volatility models has 
Leung and Sircar ( 2015l i. In practice, significantly 


fewer strikes are available for LETF options, some with wider bid-ask spreads, as compared to 
the nonleveraged counterparts. This phenomenon is partly due to the difficulty to hedge LETF 
options, which in turn impedes active market making that typically narrows bid-ask spreads. Hence, 
we consider the static hedging problem for an LETF option using options on the reference index 
(nonleveraged underlying). With the optimal static portfolio, we obtain a candidate practical 
hedging strategy for LETF options, and can also illustrate how the LETF option price and payoff 
relate to those of the vanilla options 

A call option written on the LETF L with leverage ratio i has the terminal payoff 


Et = {Lt - K')'^ = 


^Loexp [i{X^ - X^) + - K'Y =: h{K',X^, {X^)t). 


We wish to statically hedge this option using a bond, and European call and put options written 
on the reference index S. It is also possible to consider hedging an LETF option by dynamically 
trading the reference index provided it is liquidly traded, though this is beyond the scope of this 
paper. 


Hedging with discrete strikes (ATj) and a bond 

We consider the setting of Section [3^ Specifically, we assume that European-style calls/puts are 
available at a discrete number of strikes {Ki)i^i. We set Zt{0) = Bt = 1 and for i 7 ^ 0 we set 
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ZT{i) = g{Ki,ST)- For simplicity, we consider the case of no cost constraints. In order to obtain 
the optimal static portfolio vr* using Theorem 13.61 the expectations in (I3.26|] we must compute are 


= 1, V’o,i = '0i,o = E [fif(iFi,e^r)] 

= 7o = E[h(i^',X^,(Xi)T)], (5.6) 

= W.[g{Ki,e^T)h{K'{X^) t)], i,j ^ 0. 

Remark 5.1. If has constant drift and volatility ai, then 

X^ = X^ + niT + uiIFt, and {X^)t = alT. 

Also, the expectations in (|5.6p can be computed explicitly. In fact, since {X^)t = (JiT is determin¬ 
istic, we see from ()5.5I) that an option on L is simply a power option on S, which is also perfectly 
replicable with a static portfolio with a bond, a forward and calls/puts on S at all strikes K G [0, oo) 
without cost constraint. 


We wish to cons ider an i ncomp lete market model under which {X^)t is stochastic. This leads 


us to work with the 


HestonI (Il993l i model: 


dXl = {m- \X^) dt + Y^dW.S 

dAl = K{e - X^)dt + 5^^{pdWl + ^/^^IT2). 


(5.7) 


While the joint density of {X^, {X^)t) is not available, the joint characteristic moment-generating 
function of {X)p, {X^)t) is well known. We define 


^'(r,xi,X 2 ;^, A) :=Ea;i,a ;2 [exp (i^XyA(A^) t)] . 


An explicit expression for T is given in (jDrimusl . l2012l . Proposition 2.1). For a suitable function / 
we have 

OO ^?7+ioo ^ 


E,,,, J(Xi, (X1 )t) = 


1 


27ri 


de 


d7y^'(r,xi,X2;C, A)/(^, A), 


OO ^77—100 


Jn Jr+ 


(5.8) 

(5.9) 


where r/ > 0 is a positive constant chosen to the right of any singularities in the integrand of (j5.8p . 
Note, in certain cases, one may need to fix an imaginary component of ^ in (15.911 in order for the 
Fourier-Laplace transform to converge. For example, consider the function 




£,0>O. (5.10) 
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Inserting (15.101) into (15.9p and integrating yield 


The Fourier-Laplace transform / allows us to compute the Fourier-Laplace transforms of 

h{K\x,y), g{Ki,e^), h{K' ,x,y)g{Ki,e") and g{Ki,&^)g{Kj,&^) 

since these functions share the same form as / in (|5.1U|) . Inserting (15.IIP into (15.8|) and integrating, 
one can compute numerically the values of the expectations appearing in (|5.6P . 

Figure [3] illustrates the optimal static portfolio weights Ti*{Ki) for both the triple long LETF 
(£ = +3) and the triple short LETF = —3), as well as the portfolio profile i{ST) given by (15. 4p . 
If the integrated variance is very small, then we see from (15.5P that log(L'r/.ho) ~ £ log (St/S o). 
Thus, if (. is positive, then L moves in the direction of S and if £ < 0, then L moves in the direction 
opposite of S. Since the claim to be hedged Et = {Lt — K')'^ has a positive payoff if and only if 
Lt > K', we expect that, for £ > 0, the optimal hedge vr* would be to hold calls on S with strikes 
above K', and for £ < 0 the optimal hedge vr* would be to hold puts on S with strikes below K'. 
And, indeed, this is precisely what we observe in Figure [3l Moreover, we notice that the hedging 
portfolio for the (—3)-LETF call involves many long positions in puts on the reference index. This 
is not surprising given the fact that calls on the (—3)-LETF and puts on the reference index are all 
bearish positions. On the other hand, for hedging the (+3)-LETF call, the static portfolio weights 
exhibit oscillating behavior with more long positions in OTM calls. In both cases, the weights are 
assigned heavily on options that realize a positive payoff when the LETF options do also. 

5.3 Hedging a geometric Asian call with European options 

As our last example, we consider static hedging another path-dependent derivative. First, let 
5 = where 

dXl = g{t,Xl)<lt + a{t,Xl)<lWu 

That is, X^ has local volatility dynamics. We fix a time T > 0 and introduce process A = e^^ 
where satisfies 

dA/ = ^A/dt, Xl = 0. 

Note that A|. is the time average of A^ over the interval [0, T]. We consider a geometric Asian call 
with the terminal payoff: 

Ht = {At - K')+ = (e^r _ i^')+ =; h(A', A|). 
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Continuum of Strikes K G [L, R) and a bond 

Following the setting of Section [3] without the cost constraint, we construct a static portfolio 
using calls/puts available at every strike K G [L,R), plus bonds, if needed. In this example, the 
expectations that need to be computed in order to implement Theorem 13.21 are 

S = E[(e^r _e^)2], ,5 = E[(e^T _e^)], z{K) = E[g{K,e^T)]^ (5.12) 

yiK)=E[ie^T-e-)g(K,e^h], ^ = E[hiK', X^)], 9 = E[{e^T - e-)hiK', X^)], 

along with the ratio that appears in ()3.13l) . 

If the drift y and volatility a of are constant, then is jointly Gaussian. The mean 

vector and covariance matrix of {X^,X^) are 

) ’ ‘ ’ U(<r-'')2r ' 

The approximate correlation is p = \/3/2 ps 0.866. In fact, all the expectations in (j5.12p can 
be computed explicitly. Since we have already computed optimal hedges for correlated Gaussian 
random variables in Section EH we shall consider a sophisticated model next. 

We now suppose that X^ admits the Gonstant Elasticity of Variance (GEV) dynamics 

dXj = (m- l 52 g 2 (r,-i)xi^ ^ (5.13) 


.0 


The density of X^ is known in the GEV setting due to results from ICoxI (jl975l ). though the joint 
density of (V^, X^) has to be computed numerically. We follow the approximation methods outlined 
in Section E] to compute the optimal static hedges. 

In Figure m we plot the optimal unconstrained static hedging portfolio 7r*{K) as a function of 
K. Notice that the static portfolio includes most options around the strike K' of the Asian call. 
Moreover, most of the options are held long, except for the small short positions taken for the deep 
OTM calls (for K > 1.25K' in the figure). The profile resembles that in the case of hedging with 
options on a correlated asset in Figure [U (top left). This is not surprising given that S and its 
geometric average A are positively but not perfectly correlated. We also plot the hedging portfolio 
terminal value <I>(5r) as a function of St (see (I5.3p i. and it appears to be increasing convex, similar 
to a call payoff as expected. However, the slope of the payoff is significantly less than 1 even for 
large St because the Asian call payoff depends not only on St but the average of S over time. 


6 Conclusion 


We have considered the static hedging problem using a bond, a forward and a (possibly semi¬ 
infinite) strip of calls/puts, and a subset of these instruments. The optimal strategy is derived 
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wi thout assuming dynami cs of the underlying asset(s), and is shown to reconcile with the result 


of 


Carr and MadanI (Il998l ). which is a special case of our framework. A useful connection is estab¬ 


lished between the optimal static hedging strategy with discrete option strikes and our proposed 
continuous-strike strategy implemented over discrete strikes, and we show the contrasting behav¬ 
iors of these strategies. The numerical implementation of the optimal strategy is conducted in a 
general incomplete Markov diffusion market, with examples of exotic/path-dependent claims such 
as options on a non-traded asset, as well as Asian and LETT options. 

For future research, one direction is to generlize our results for static hedges up to a random time, 
rather than a fixed terminal time considered here. This extension will provide new alternative ways 
to statically hedge claims such as barrier options and American options. I n other related direct ions, 
one can consider semi-static hedges with finite number of rebalancing (see Carr and wJ (2013)), or 


incorporate static positio n s into portfolio optimization problems; see, for example, 
( 2005 1: Leung and Sircar (2009). 


Ilhan and Sircar 
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A Proof of Proposition 13.6 

Optimization problem (|3.25l) is a finite-dimensional quadratic programming problem with linear 
constraints. By Assumption 13.51 the matrix -0 is positive definite, and thus, invertible. Define the 
Lagrangian 


L(7r, A) := J(7r) - A • {H{Tr) - C). 


The Karush-Kuh n-Tucker (KKT) conditions, w hich are necessary and sufficient for this convex 
optimization (see IBovd and Vandenberghel (j2004l lL are 


stationarity : 


complementary slackness : 


0 = d^^L{7T,\) 

i j i i 

= 2 ^ - 2 , V/s, (A.l) 

0 = A • (iL(7r) - C) 
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In matrix notation, equation (jA.ll) becomes 


V’TT = 7 + —Z 


vr = 7r(A) ='iIj ^ ^7 + (A.3) 

TTi = TTiiX) = (^7J + . (A.4) 


Inserting expression (IA.4I) for vrj into (IA.2P we obtain 


0 = A- 




A.(5T^-i + 


Equation (|A.5ll has two solutions: 


A = 0, 


and 


Inserting expression (jA.6l) into (IA.3I) yields (I3.27p . 


fCjzlYLh'] 


(A.5) 


(A.6) 
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Figure 1: Left: For unconstrained static hedging problem, we plot the optimal density Tr*{K) is plotted 
over K of puts/calls used in the unconstrained static hedging problem (top left), and the terminal portfolio 
profile ^{St) according to (15.31) (bottom left). The dashed, dot-dashed, and dotted lines correspond to 
p = (0.9,0.7,0.5), respectively. The solid line represents the payoff to be hedged, {Vt — K')^. The parameters 
used are pi = p .2 = 0.1, ai = a 2 = 0.2, So = Vq = K' = 1 and T = 0.5. Right: Let c be the cost of the 
unconstrained optimal portfolio. We plot the optimal density •k*{K) (top right), and the portfolio profile 
<i>(5'T), given by (15.3|) . both for the cost-constrained static hedging problem. The dashed, dot-dashed, and 
dotted lines correspond to cost constraints of C = (c, 0.75c, 0.5c), respectively. On the right panel, the 
parameters used are pi = P 2 = 0.1, tri = CT 2 = 0.2, So = Vo = K' = 1, p — 0.55 and T = 0.5. 
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Figure 2: Left: For the unconstrained static hedging problem with discrete hedging assets, we plot the 
optimal number TT*{Ki) of puts/calls over discrete strikes. The black, gray and white bars correspond to 
p = (0.9, 0.7,0.5), respectively (top left). The portfolio profile $(S't), given by (15.4p . is plotted as a function 
of St (bottom left). The dashed, dot-dashed, and dotted lines correspond to p = (0.9, 0.7, 0.5), respectively. 
The solid line corresponds to the payoff to be hedged, (Vr — 77')+, plotted as a function of Vt- For the 
top/bottom left panels, the default parameters are pi = P 2 = 0.1, cti = (72 = 0.2, = Vq = 1 and 

T = 1.0. Right: Let c be the cost of the unconstrained optimal portfolio. We plot the optimal density 
TT*{Ki) of puts/calls for the cost-constrained static hedging portfolio, corresponding to cost constraints 
C = (c, 0.75c, 0.5c) (black, gray, and white bars, respectively) (top right). The terminal value $(S't) for the 
cost-constrained optimal portfolio as in (15.41) is shown for cost constraints of C = (c, 0.75c, 0.5c) (bottom 
right). For the top/bottom right panels, the default parameters are pi = P 2 = 0.1, ai = a 2 = 0.2, 
So = Vo = I, p = 0.9 and T = 1.0. 
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Figure 3: In the Heston model (15.7|) . we consider statically hedging a triple long LETF (left, £ = +3) and 
triple short LETF (right, i = —3) by including bonds, and calls/puts on the reference index at discrete strikes 
{Ki)i^j. Top: We plot as a function of Ki the optimal number Tr*{Ki) of puts/calls in the unconstrained 
optimal static hedging portfolio. Bottom: The portfolio profile $(S't) as a function of St (see (15.41) 1 for the 
two cases. The parameters used are xi = 0, X 2 = 0.04, Lq = K’ = I, m = 0.1, 0 = 0.04, k = 1, p = 0, 
5 = 0.1 and T = 0.25. 
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Figure 4: In the CEV model (15.131) . we consider statically hedging an Asian call using bonds, forwards, 
and calls/puts on the underlying stock at all strikes K G [0, oo). For the put and calls in the unconstrained 
static hedging portfolio, we plot the optimal density •k*{K) as a function of strike K (left). Also, we plot the 
portfolio profile 4 >(S't) as a function of St (right), according to (15.31) . The parameters used here are xi = 0, 
K’ = l,e = 0.04, m = 0.1, 5 = 0.2, rj = 0.7, and T = 1.0. 
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